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Abstract 

Recently we have derived an effective one-dimensional nonpolynomial 
Schrodinger equation (NPSE) [Phys. Rev. A 65, 043614 (2002)] that accu- 
rately describes atomic Bose-Einstein condensates under transverse harmonic 
. confinement. In this paper we analyze the stability of the solitary-wave so- 

o . 

lutions of NPSE by means of the Vakhitov-Kolokolov criterion. Stable self- 

O 

£^ . focusing solutions of NPSE are precisely the 3D bright solitons experimentally 



> 



found in Bose-Einstein condensates [Nature 417, 150 (2002); Science 296, 
1290 (2002)]. These self-focusing solutions generalize the ones obtained with 



the standard nonlinear Schrodinger equation (NLSE) and take into account 
the formation of instabilities in the transverse direction. We prove that nei- 
ther cubic nor cubic-quintic approximations of NPSE are able to reproduce 
the unstable branch of the solitary-wave solutions. Moreover, we show that 
the analitically found critical point of NPSE is in remarkable good agreement 
with recent numerical computations of 3D GPE. We also discuss the formation 
of multiple solitons by a sudden change of the sign of the nonlinear strength. 
This formation can be interpreted in terms of the modulational instability 
of the time-dependent macroscopic wave function of the Bose condensate. In 



particular, we derive an accurate analytical formula for the number of multiple 
bright solitons. 
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I. INTRODUCTION 



A nonlinear Schrodinger equation with nonpolynomial nonlinearity, the so-called non- 
polynomial Schrodinger equation (NPSE), has been recently introduced and studied [1-3]. 
This equation has been derived in the context of Bose-Einstein condensation [4,5] from the 
3D Gross-Pitaevskii equation [6,7] as a reliable effective ID equation which describes the 
axial wavefunction ip(x,t) of a cigar-shaped condensate confined in the transverse direction 
(y, z) by a harmonic trapping potential. NPSE has been successfully used to numerically 
investigate [8] the dynamics of Bose condensate solitary-waves experimentally observed by 
two groups [9,10]. In this letter we rigorously determine the linear stability domain of 
the solitary-wave solutions of NPSE by using the Vakhitov-Kolokolov criterion [11,12]. In 
particular, we find that NPSE admits both stable and unstable stationary solitary-wave 
solutions. 



II. GENERALIZED SCHRODINGER EQUATIONS 



For the generalized nonlinear Schrodinger equation (GNLSE), given by 

8 1 d 2 



1~^r~ + 



+ F 



^ = 0, 



(1) 



dt 2 dx 2 

where F{\ip\ 2 ) is the nonlinear term, the linear stability of solitary- waves is determined by 
the Vakhitov-Kolokolov criterion [11-13]. With the standard position 



1>{x,t) = fa) e**, 



(2) 



where uj is the frequency of the fase and <j){x) is a real scalar field, from GNLSE one obtains 
the stationary second-order differential equation 

"1 d 2 



2dx 2 



+ F{4> 2 ) 



(j) = uj <p . 



(3) 



The first derivative N' s {uj) of the number of particles invariant 

N s (uj) = j dx 4> 2 {x) , 
3 



(4) 



calculated for the soliton solution <f>(x), allows one to predict the parameter region in uj 
where the soliton amplitude can grow or decay exponentially with a nonzero growth rate. 
The solitons are stable if N' s {uj) > and unstable if N' s {uj) < [13]. 



In the case of the self-focusing nonpolynomial Schrodiger equation (NPSE) [1-3] the 
nonlinearity of Eq. (1) and Eq. (3) is given by 



where g is the strength of the attractive interaction between particles in the Bose-Einstein 
condensate. In particular, the parameter g is related to the s-wave scattering length a s 
and to the characteristic length a±_ of the transverse harmonic confinement by g = 2\a s \/a± 



Given a GNLSE like Eq. (3) it is quite natural to try a power expansion of F(<p 2 ). In 
the case of NPSE one finds 



The lowest order in g<f) 2 gives the familiar one-dimensional self-focusing cubic NLSE. By 
including the next term of the power expansion one has a cubic-quintic NLSE, which has 
been investigated in several papers (see for instance [13]). But, as we shall show, the cubic- 
quintic approximation of NPSE is not able to correctly describe the stability of NPSE solitary 
waves. 

The exact stationary solitary-wave solutions of NPSE can be obtained in the following 
way. A simple constant of motion of the stationary NPSE is 



III. NONPOLYNOMIAL SCHRODINGER EQUATION 





[1-3,8]. 



F(0 2 ) = -l+^ 2 + ^ 2 4 + .... 



(6) 




(7) 



where 



4 



(8) 



is the effective potential energy of the system. In Figure 1 we plot the potential energy V{4>) 
of the full NPSE and its cubic and cubic-quintic approximations. As shown by Figure 1, 
for uj < V{4>) is a double- well potential and the two approximations are quite close to the 
exact curve, apart near = ±l/yfg, the end-points of the NPSE which are related to the 
the singularity of the term (1 — gcj) 2 ) 1 / 2 . By imposing the boundary condition 0^0 for 
x — > oo one finds E = and so 



The previous formula gives the integral equation 



(9) 



J dx = j d(f) 



(10) 



from which one obtains the solitary- wave solution <f>(x) written in implicit form (see also [8]) 
1 



V2x = 



arcth 



VI - g<j)' 2 + oo 

1+LU 
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Vl - g<t> 2 + u 

l-UJ 



(11) 



IV. VAKHITOV-KOLOKOLOV STABILITY CRITERION 



The number of particles N s (uj) of Eq. (4) can be calculated by observing that 



dx 



from which one gets 



N 3 (u) = / # 







y/2y/l - get) 2 + 2lu 
This integral is analytically solvable and one obtains 

2^2 



N s (u) 



3^ 



(1 - 2u) Vl +u) . 



(12) 



(13) 



(14) 
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In Figure 2 we plot the function N s (u) for some values of the nonlinear strength g. The 
Vakhitov-Kolokolov criterion [11,12] is based on the study of the sign of the first derivative 

iV» = (18) 
g vl +co 

of the number of particle N s (u). One easly finds that the solitary- wave 4>(x) given by Eq. 
(11) is linearly stable for — 1 < uo < —1/2 and unstable for —1/2 < uo < 1/2. Moreover 
N s (uo) has its maximum value N™ ax = 4/(3g) at the critical frequency uo = —1/2. 

In Figure 2 we also plot N 3 (u) obtained from the cubic NLSE and the cubic-quintic 
NLSE. In the case of cubic NLSE one finds 

N a (u) = ^y/T+Z, (16) 

from which one immediately obtains that the solitary-wave solutions are always stable. In 
the case of the cubic-quintic NLSE the number of particles is instead given by 



2 f TV 

N s (u) = - - - arctg 



V2 



2^/T+ 



UO 



(17) 



9 V 2 

Remarkably also in this case the Vakhitov-Kolokolov stability criterion tells us that there 
are no unstable solitary-waves. Thus the quintic term in the power expansion of NPSE is 
not sufficient to determine the stability domain of NPSE. 

As previously stated, in the context of Bose-Einstein condensation [1-3] the nonlinear 
strength g is given by g = 2|a s |/a_i_. By using this expression one finds that the maximum 
value of N s (lu) for the full NPSE is given by 

N™ ax = -^r. (18) 
3|a s | 

This value is very close to the critical number N c of attractive Bose-condensed atoms under 
transverse harmonic confinement found in [14] with a full numerical solution of the 3D GPE: 

N c = 0.676^ . (19) 
\a s \ 

Beyond iV c there is the so-called collapse of the Bose-Einstein condensate. The collapse is 
due to the fact that the solitary- wave solution becomes linearly unstable, i.e. N s (u) < 

6 



in our NPSE. Note that this collapse has nothing to do with the singularity of the term 
(1 — gcf 1 ) 1 / 2 that is present in NPSE: in fact, also at the collapse the condition g<j) 2 < 1 is 
satisfied. 



V. MODULATION AL INSTABILITY 



In a recent experiment [9] it has been reported the formation of bright solitons in a Bose- 
Einstein condensate of 7 Li atoms induced by a sudden change in the sign of the scattering 
length from positive to negative. The formation of these solitons can be explained as due to 
the modulational instability of the time-dependent wave function of the Bose condensate, 
driven by imaginary Bogoliubov excitations [15]. 

Let us consider a Bose condensate with a repulsive inter-atomic interaction (a s > 0). 
The condensate is described by Eq. (1) and Eq. (5), setting g' = —g instead of g in Eq. (5). 
Under box axial confinement, the stationary homogeneous wave function is given by 



for —L/2 < x < L/2 and zero elsewhere. L is the size of the confining box and N the total 
number of bosons. The Bogoliubov elementary excitations e k of the static Bose condensate 
4>(x) are found by looking for solutions of the form 



and keeping terms linear in the complex functions u(x) and v(x) (linear-stability analysis). 
In the quasi-lD limit (defocusing NLSE) one finds 



where ui = —g'n, n = N/L is the axial density and g' = 2a s /a±. 

By suddenly changing the scattering length a s to a negative value, the excitations fre- 
quencies corresponding to 




(20) 




(21) 




(22) 



7 



become imaginary and, as a result, small perturbations grow exponentially in time. This 
penomenon is known as modulational instability [16,17]. It is easy to find that the maximum 
rate of growth is at k = k c /y/2. The wavelength of this mode is Ao = 27r//c and the ratio 
L/Aq gives an estimate of the number Nbs of bright solitons which are generated: 



The predicted number Nbs of bright solitons is in very good agreement with the numerical 
results of 3D GPE [18] and in rough agreement with the experimental results [9]. Note that 
the finite resolution of the imaging process reduces the number of detected solitons [18]. 

Finally, we observe that Eq. (24) has been obtained by using NLSE. NPSE gives that 
same results of NLSE in the quasi-lD limit (of experiment [9]) but predicts the formation 
a large number of solitons also close to the critical point gn — 1. The consequences of this 
prediction will be discussed elsewhere. 



In this paper we have shown that the solitary- waves studied in [8] belong to the stable 
branch — 1 < uo < — 1/2 we have here analytically found by using the Vakhitov-Kolokolov 
criterion. In this branch, growing uj the number of particles N s increases and for a fixed N s 
the solitary-wave solution is fully determined. Since the study of the 3D model is usually 
time consuming, it is nice to have a simple ID model to understand the dynamics of a Bose- 
Einstein condensate. NPSE has the same numerical complexity of other one-dimensional 
nonlinear Schrodinger equations but we have shown that it is not well approximated by the 
cubic or cubic-quintic Schrodinger equations. NPSE reflects correctly the 3D dynamics of 
Bose condensed solitary-waves and it predicts with considerable accuracy the collapse of the 
condensate. Moreover, we have shown the NPSE explains remarkably well the formation 
of multiple bright solitons by a sudden change of the sign of the nonlinear strength from 




(24) 
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CONCLUSIONS 



positive to negative, as a consequence the modulational instability of the time- dependent 
wave function. 
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FIG. 1. Effective potential energy V{<j>) of the NPSE. Exact: full line; cubic approximation: 
dotted line, cubic-quintic approximation: dashed line, uj is the phase frequency and g is the 
nonlinear strength. 
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FIG. 2. Number of particles N s (uj) as a function of the phase frequency uj for the NPSE. 
Exact: full line; cubic approximation: dotted line; cubic-quintic approximation: dashed line, g is 
the nonlinear strength. 
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